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A new quantum algorithm is proposed to solve Satisfiability(SAT) problems by taking advantage 
of non-unitary transformation in ground state quantum computer. The energy gap scale of the 
ground state quantum computer is analyzed for 3-bit Exact Cover problems. The time cost of this 
algorithm on general SAT problems is discussed. 
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Quantum computer has been expected to outperform 
its classical counterpart in some classically difficult prob- 
lems. For example, the well-known Shor's factoring 
algorithm^ and Grover's algorithm^ accelerate expo- 
nentially and quadratically compared with classical al- 
gorithms. It is a challenge to find whether quantum 
computer outperforms its classical counterpart on other 
classically intractable problems 0, Q, which cannot be 
solved classically in polynomial time of N, the number 
of input bits. Especially interesting are the NP-complete 
problems Q, which include thousands of problems, such 
as the Traveling Salesman problem^ and some satisfi- 
ability (SAT) problems. All NP-complete problems can 
be transformed into each other by polynomial steps. If 
one of the NP-complete problems can be solved in poly- 
nomial time by an algorithm in the worst case, then all 
NP-complete problems can be solved in polynomial time. 
However, it is widely believed that such a classical algo- 
rithm doesn't exist. 

In this paper we explore the idea of ground state quan- 
tum computer (GSQC) 0, H Hi H(i| . and propose a new 
algorithm to solve SAT problems. A iC-SAT problem 
deals with N binary variables submitted to M clauses 
with each clause Ci involving K bits, and the task is to 
find TV- bit states satisfying all clauses. When K > 3, K- 
SAT is NP-Complete, and some instances become clas- 
sically intractable when the parameter a = MJN , as 
M, N -> oo, is close to threshold a c (K)\L& 

A standard computer is characterized by time depen- 
dent state as: \^{U)) = Ui\ij)(ti-l))i where U denotes 
instance of the i-th step, and Ui represents for unitary 
transformation. For GSQC, the time sequence is mim- 
icked by the space distribution of the ground state wave- 
function \tpo). 

As proposed by Mizel et.al.0, a single qubit may 
be a column of quantum dots with multiple rows, and 
each row contains a pair of quantum dots. State |0) 
or |1) is represented by finding electron in one of the 
two dots. GSQC is made up by circuit of multiple 
interacting qubits, whose ground state is determined 
by the summation of single qubit unitary transforma- 
tion Hamiltonian h(Uj), two-qubit interacting Hamilto- 
nian h(CNOT), boost Hamiltonian h(B,X) and projec- 
tion Hamiltonian h(\"f), A). The energy gap between the 
ground state and the first excited state determines the 
efficiency of GSQC[]3. 



The single qubit unitary transformation Hamiltonian 
has the form: 



h j {Uj 



C}_ x Cj- X + C}C. 



UjCj-i + h.c. 



,(1) 



where e defines the energy scale of all Hamiltonians, 

Cj = Cj Cj x , cj is the electron creation operator 

on row j at position 0, and Uj is two dimension matrix 
representing for unitary transformation from row j — 1 
to row j. The boost Hamiltonian is: 



h j {B, A) = e 



C 



h.c. 



,(2) 



which amplifies the jth row wavefunction amplitude by 
large number A compared with (j — l)th row in |^o) • The 
projection Hamiltonian is 



ft* (| 7 >,A) = e 



1 

-1,7 + 



\ ( c L c i-^ + h - c ) 



where \j) represents for state to be projected to on row j 
and to be amplified by A. The interaction between qubit 
a and f3 can be represented by h(CNOT): 

h j a ^CNOT) = eOlj^C^-x^Cpj 

+K( I ) c l, j -i G 0,j-i + c L-,o c <*,j,o W 
+4j,i^J,ifcj(tf)- ( 4 ) 

where for c\ b c , its subscription a represents for qubit a, b 
for the number of row, c for the state \c). All above men- 
tioned Hamiltonians arepositive semidefinite, and are 
the same as those in 0, la 13 • Only pairwise interaction 
is considered for interacting Hamiltonians. 

The input states are determined by the boundary con- 
ditions applied upon the first rows of all qubits, h — 
E(I + J2i a i a i)i with Oi being Pauli matrix, J2i a i = 1 
and E being large compared with e|lfj|. 

To implement any algorithm, on final row of each qubit 
boost or projection Hamiltonian is applied so that |?/>o) 
concentrates on the position corresponding to the final 
step in standard paradigm, hence measurement on GSQC 
can read out desired information with appreciable prob- 
ability. 
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As shown in [13 GSQC circuit may have exponentially 
small energy gap depending on detail of circuit, and as- 
sembling GSQC circuit directly following algorithm for 
standard paradigm, such as quantum Fourier transform, 
leads to exponentially small energy gap. In order to avoid 
small gap, teleportation box, as shown in Fig.QJ, is in- 
troduced on each qubit between two control Hamiltoni- 
ans. The teleportation boxes make all qubits short (the 
longest qubit has length 8), on the other hand, for ar- 
bitrary GSQC circuit they make the energy gap only 
polynomially small A cx e/A 8 if all boost and projec- 
tion Hamiltonians have the same A value. To determine 
magnitude of A, one only needs to count the total num- 
ber of qubits in the circuit, which is proportional to the 
number of control operation in an algorithm, say N k , 
then the probability of finding all electrons on final rows 
is P ~ (1 — C/X 2 ) N with C being 8, the maximum 
length of qubit. In order to have appreciable P, we set 
A « iV fc / 2 , hence A cx e/N Ak . The details can be found 
in 0. 

While a time-dependent standard quantum computer 
makes unitary transformation from one instance to the 
next, GSQC may have non-unitary transformation from 
one row to the next, such as the boost Hamiltonian 
h(B, A) and projection Hamiltonian h(\j), A). Especially 
the projection Hamiltonian, which mimics measurement 
in standard paradigm, can amplify the probability of cer- 
tain state to be "measured", hence GSQC owns advan- 
tage over standard quantum computer. 

A simple example, although of no practical interest, 
demonstrates this advantage: to teleport quantum state 
from qubit 1 to qubit 2, then to qubit 3, and so on to 
qubit N. By standard quantum computer, the probabil- 
ity to successfully realize this series of teleportations is 
(1/4) N because each teleportation process only has 1/4 
probability to succeed [l4|: while by GSQC, the probabil- 
ity is (1 - 8/\ 2 ) 3N : setting A = s/T)N, then P » e - 2 ^ D , 
and energy gap is A cx e/(D 4 N A ) 10}. Thus GSQC only 
costs polynomially long time to finish the task while stan- 
dard paradigm needs exponentially long time. 

The advantage of GSQC makes new quantum algo- 
rithm possible. Here I present a quantum algorithm to 
solve SAT problems as shown in Fig.©, a GSQC circuit 
to solve a 3-SAT problem with only 9 bits. It's easy to 
extend to A-bit problem. Each clause is implemented by 
a "filter box" , and the circuit inside each filter box makes 
sure that on rows immediately below it the states satisfy- 
ing clause Ci have much larger amplitudes than other un- 
satisfying states, or we can say those unsatisfying states 
are filtered out. This can be realized by projection and 
boost Hamiltonians, and the detail will be given in the 
following example. In the figure, the input state on the 
top row is (|0) + |1»(|0) + |1»..-(|0) + |1», which is de- 
termined by the boundary Hamiltonian, h = E(I — a x ); 
the clause involving qubit 1, 2 and 3 is implemented by 
filter box 1, the clause involving qubit 2, 3 and 4 imple- 
mented by filter box 4, the clause involving qubit 3, 4 
and 8 implemented by filter box 6, etc. 



P(0) |P(0) 

(!) 



FIG. 1: Design for "teleportation box". The circuit inside 
box is similar to the teleportation circuit in Q. Label / rep- 
resents for identity transformation Hamiltonian h(I), H for 
Hadamard transformation Hamiltonian h(H) and P(0) for 
projection Hamiltonian h(\0),X). 
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FIG. 2: Design for circuit solving SAT problem with each 
clause involving 3 bits. Box labeled "Filter Box" represents 
for circuit as shown in Fig, d^l . 



When all constraints are implemented, at ground state 
the states measured on the final rows of the N qubits 
should be superposition of states satisfying all con- 
straints. I will show no backtracking is needed later. 

Now I give an example on how to implement a filter 
box. We focus on the 3- bit Exact Cover problem^, an 
instance of SAT problem, and belongs to NP-complete. 
Following is definition of 3-bit Exact Cover problem: 

There are N bits zi, Z2, Zff, each taking the value 
or 1. With O(N) clauses applied to them, each clause is 
a constraint involving three bits: one bit has value 1 while 
the other two have value 0. The task is to determine the 
N-bit state satisfying all the clauses. 

The algorithm is implemented by the circuit in Fig.J2J). 
Considering any one of the clauses, in GSQC a filter box, 
involving three qubits i, j and k, which are represented 
by gray dot columns in Fig.ljJJ, we add two ancilla qubits: 
qubit 1 and qubit 2, which are represented by dark dot 
columns. Qubits i, j, k at the first row are in the state 
(|1) + |0)) if they have not experienced any clause yet, 
and the two ancilla qubits are in state |0) and |0) on top 
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FIG. 3: Design for "filter box". The labels on the lines repre- 
sent for corresponding Hamiltonians: / for h(I), CNOT for 
h(CNOT), P(l) for projection h(\l},\) et. al. At the final 
rows, B/I represents for boost Hamiltonian h(B, A) if there is 
no more clause to be applied to this qubit, otherwise, repre- 
sents for identity Hamiltonian h(I). There are teleportation 
boxes, not shown in figure, inserted on all qubits between two 
control Hamiltonians. Some dots marked by X or Y are for 
demonstration convenience in text. 



rows by boundary Hamiltonians, where \j) corresponds 
to state of ancilla qubit 1, and I7) to state of ancilla qubit 
2. 

Inside the dashed triangle of Fig.©, after the first 
CNOT, we obtain state |i)|l) + |6)|0); after the second 
CNOT: |1)|1)|0) + |6)|0)|0) + |6)|1)|1) + |1)|0)|1); after 
the third CNOT: 



|i)(|i)|o)|o) 
|6)(|i)|i)|o) 



|0)|1)|0) 

|o)|i)|i) 



|0)|0)|1) 

|i)|o)|i) 



|0)|0)|0)). 



Immediately below the triangle, if the system stays at 
ground state, electron in ancilla qubit 1 is measured to 
be on the row labeled by X and at state |1), and the three 
electrons on qubit i, j, k are all found on the rows labeled 
by X, then the three-qubit states satisfy the clause except 
for |1)|1)|1>. 

The ancilla qubit 2, starting at state |0), experiences 
CNOT gates controlled by qubits j and k, and J2(±7r/4) 



is defined in [3 as R y {±n/4), as shown within the dotted 
pentagon in Fig. © . All those transformations happened 
inside the dotted pentagon are equivalent to a Toffoli gate 
except for some unimportant phases [l5j: if both qubits 
j and k are in state |1), then the ancilla qubit 2 reverses 
to state |1), otherwise, it remains at state |0). After this 
nearly Toffoli transformation, if at ground state electrons 
in qubit j, k and ancilla qubit 2 are found on rows labeled 
by Y, and ancilla qubit 2 is at |0), then the three qubits 
will be at |0)(|0)|0) + |1)|0) + |0)|1)). Thus if at ground 
state all electrons are found on rows immediately below 
both the dashed triangle and the dotted pentagon, and 
if ancilla qubit 1 is at |1) and ancilla qubit 2 at |0), then 
the three qubits i, j, k and two ancilla qubits will satisfy 
the clause: 



|1)|0)(|1)|0)|0) + |0)|1)|0) + |0)|0)|1)). 



(5) 



In order to make the right states pass through the filter 
box with large probability, we add projection Hamiltoni- 
ans and boost Hamiltonians as shown in the lower part 
of Fig.©. The projection Hamiltonians on final rows of 
two ancilla qubits limit and amplify them at the states 
we prefer: ancilla qubit 1 at |1), and ancilla qubit 2 at |0). 
If a qubit doesn't experience any more clause, it will end 
with a boost Hamiltonian, otherwise, its quantum state 
will be teleported to a new qubit through teleportation 
box, not shown in Fig.©, and the new qubit continues 
experience more clauses. Thus the projection Hamilto- 
nian on two ancilla qubits and boost Hamiltonian on the 
three qubits make sure that the ground state wavefunc- 
tion concentrated on the final rows in Fig.© with state 
at Eq.©. 

Noting that in the filter box all the three qubits i, j, 
and k always act as control qubits, thus the entanglement 
of these three qubits with other qubits not involved in 
this particular clause still keep the same. When adding a 
clause, the resulted states satisfying this clause will also 
satisfy all previous applied clauses. Thus unlike classical 
algorithm, no backtracking is needed. 

In the circuit of Fig.©, if there is at least one solution, 
and all electrons are simultaneously found on the final 
rows of all qubits, then the reading of the JV-bit states 
satisfying all constraints. 

In order to keep the energy gap from being too small, 
like in |lOj . on every qubit teleportation boxes are in- 
serted between two control Hamiltonians, thus the to- 
tal number of qubits increases while the energy gap 
A cx e/A 8 if in all boost and projection Hamiltonian the 
amplifying factors have the same value A. 

For one clause, or a filter box, it needs 10 teleporta- 
tion boxes (each teleportation box adds two more qubits) 
on the original five-qubit circuit, noting that on the end 
of qubit i, j and k in Fig.© teleportation boxes are 
needed because more clause will be added. Thus adding 
one more filter box means adding 20 more qubits. The 
number of clause for a NP hard 3-bit Exact Cover prob- 
lem is about the same order as the number of bits N 
say aN with a being 0(1), then there are totally 20aN 
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qubits and each of them ends with either projection or 
boost Hamiltonian. Probability of finding all electrons 
at the final rows is approximately 



P»(l- C/A 2 ) 



20aN 



(6) 



where C — 8, the length of the longest aubit[Tol|. It is as- 
sumed that, at ground state, in each filter box the ancilla 
qubit 1 and 2 have appreciable probability in |1) and |0) 
states respectively before projection Hamiltonians. We 
will address situation when the assumption is violated. 

In order to make the probability independent of num- 
ber of bits N, we take A 2 = DN, where D is an arbitrary 
number. Then as ./V becomes large, we obtain 



P « (1 - C/(DN)) 



20aN 



,-20aC/D 



and energy gapjlCJ 

Acxe/A 8 oce/(D 4 iV 4 ), 



(7) 



(8) 



from which one can estimate time cost. 

To make the GSQC circuit at ground state, we can 
use adiabatic approach: first we set A = 1 for boost and 
projection Hamiltonian on final rows of all qubits, and 
replace the single qubit Hamiltonian between the first 
two rows of all qubits by a boost Hamiltonian 



h'(B,X r ) = e 



j^cld + c\c 2 i (c\o 2 



,(9) 



so that the wavefunction amplitude of the first row is 
boosted as A' > 1. Now in the ground state the elec- 
trons concentrate at the first rows as 1/A' — > 0, thus 
the ground state is easy to be prepared, and energy gap 
A cx e/n 2 with n = 8 being the length of longest qubit. 
The next step is turning the quantity 1/A' to 1 adia- 
batically, during which the energy gap remains at e/n 2 
and the ground state wavefunction spreads to other rows 
from the first row. The third step is turning 1/A from 
1 to 1/y/DN adiabatically. In this process the energy 
gap decreases monotonically from e/n 2 to what we ob- 
tained above: e/D 4 N 4 , and wavefunction concentrated 
on the final rows of all qubit as we wish. Thus the scale 
of time cost is about T cx 1/ A 2 cx TV 8 [l^ , local adiabatic 
approach may reduce the time cost further \Tf\. 

Above analysis is under the assumption that the num- 
ber of satisfying states gradually decreases as the clauses 
are implemented one by one. There is a situation that 
might hurt our algorithm: after adding one more clause, 
if the number of satisfying states drops dramatically, our 
algorithm will be hurt. For example, if one constructs 
GSQC for the Grover's search problem with one condi- 
tion to find a unique satisfying state from 2 N states, then 



he will find that there is an ancilla qubit containing such 
unnormalized state 



|0) | satisfying) 



2 N -1 

£ 

i=l 



|1) | unsatisfying 



(10) 



before the projection Hamiltonian h(\0), A). In order to 
amplify the amplitude of the correct state on the final 
row, it requires A > 2 N / 2 , which makes the energy gap 
exponentially small. 

Does this happen to general SAT problems? In [l8| it 
was suggested that close to threshold a c computational 
complexity might be related with the forming of back- 
bone, each of a subset of bits has average value close to 1 
or in the subspace of satisfying states. The existence of 
backbone means that most satisfying states contain the 
state represented by backbone, and if adding one more 
clause kicks out the states consistent with backbone from 
satisfying subspace, the number of satisfying states drops 
dramatically. With advantage over classical algorithm, 
performance of our algorithm is not affected by forming 
of backbone, however, as more clauses applied, the dis- 
appearance of already existed backbone in the satisfying 
subspace might hurt. There is a criterion determining 
efficiency of our algorithm: the ratio Sj/Sj+i, with Sj 
being the number of solutions when the jth clause is ap- 
plied, and Sj+i the number of solutions when the (j'+l)th 
clause is applied. For example, So/ Si =8/3 for 3-bit Ex- 
act Cover problem. If Sj/Sj+i 3> 1, on the ancilla qubit 
of the (j + l)th filter box, the probability of finding elec- 
tron on its final row will be p sw (1 — CSj/(\ 2 Sj+i)). 
In order to have appreciable probability as Eq.J7|), it re- 
quires A 2 increase from DN to DNSj/Sj+i, hence the 
energy gap is also suppressed. In advance one cannot 
know what value Sj/Sj+i is, thus a overhaul factor for 
A is needed. If this ratio happens to be exponentially 
large, then our algorithm cannot solve the SAT problem 
in polynomial time. However, one might be able to iden- 
tify backbone by trials, and then choose proper order to 
implement clauses so that Sj/Sj+i can be kept small. 

In conclusion, we have demonstrated that a ground 
state quantum computer can solve a general SAT prob- 
lem. A specific example, the 3-bit Exact Cover problem, 
is given. We show that an 3-bit Exact Cover problem can 
be solved by the quantum algorithm described here, and 
the time cost is related with the number of bits N and 
parameter Sj/ Sj+\. If Sj/Sj+i stays small or only poly- 
nomially large, then the presented algorithm can solve 
this SAT problem in polynomial time. 
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